Upper bound for the conductivity of nanotube networks 
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Films composed of nanotube networks have their conductivities regulated by the junction resis- 
tances formed between tubes. Conductivity values are enhanced by lower junction resistances but 
should reach a maximum that is limited by the network morphology. By considering ideal ballistic- 
like contacts between nanotubes we use the Kubo formalism to calculate the upper bound for the 
conductivity of such films and show how it depends on the nanotube concentration as well as on 
their aspect ratio. Highest measured conductivities reported so far are approaching this limiting 
value, suggesting that further progress lies with nanowires other than nanotubes. 



The search for thin fihiis that are flexible, transparent 
and conductive is driven by their potential as transpar- 
ent electrodes [H 0, One common route is to use 
films made hy a disordered network of carbon nanotubes 
(NT) [1, S la 0, B . In this case, electrons move across 
the entire film by moving between NT in close proxim- 
ity. The conductivity is limited by the tunneling between 
tubes, which introduces a significant inter-NT junction 
resistance. To make the films more conductive one needs 
to improve the coupling between NT^which recently has 
been achieved with acid treatments [l^, |Tl| . Further at- 
tempts are being made to lower the junction resistance 
and surpass the best conductivity reported so far, which 
currently stands at cr w 6 x 10^ S/m 0, [10, [HI. 

In addition to the junction resistance, the network 
morphology also plays a role in limiting the film conduc- 
tivity. In fact, we have recently demonstrated how sen- 
sitive to the network connectivity the conductivity can 
be [T^ . This means that no matter how much progress 
is made in lowering the junction resistance, there should 
be a maximum value for the film conductivity, which is 
regulated by the network itself. This is the goal of the 
present manuscript, i.e., to obtain an upper bound for 
the conductivity of disordered NT-networks. The knowl- 
edge of this upper bound should avoid overoptimistic ex- 
pectations for the transport properties of the films. Fur- 
thermore, understanding the interplay between network 
morphology and the intrinsic conductances of NT may 
be explored to deal with films made of other nanowires 
(NW) [11,111. 

Since we are interested in the best-case scenario in 
which the electronic conductivity is at its maximum, we 
must eliminate potential sources of scattering and deco- 
herence such as structural imperfections, impurities and 
interaction with other quasi-particles. In this situation 
it is appropriate to consider a purely ballistic regime of 
transport within the NW, which calls for a quantum de- 
scription of the conductivity. NT are known to behave as 
ballistic conductors with two quanta of conductance fisj 
and are often referred to as possessing two conducting 
channels. How much interference there is between these 
channels is what determines how the network affects the 



film conductivity. Furthermore, because other quantum 
wires are populated by a different number of conduct- 
ing channels, the results here obtained are not exclusive 
for NT and should be also applicable to other materials. 
With this generality in mind, we introduce a model that 
describes the transport properties of a network of quan- 
tum wires, each one of them capable of carrying M chan- 
nels of conductance. The model consists of two parts, one 
macroscopic in which geometrical constraints define the 
connectivity of the network, and another of microscopic 
origin that accounts for the electronic structure of the 
wires. 

We start with the macroscopic part of the model by 
assuming that a film in contact to two electrodes is rep- 
resented by a number of rods of length £ and diameter 
d randomly distributed inside a rectangular box of di- 
mensions L X L X 2L. The opposing square faces of the 
box are taken as the electrodes. It is a simple task to 
count the number of rods (a^;) crossing the electrode 
walls as well as (a), here defined as the average number 
of contacts (per rod) between neighboring rods. Geo- 
metrical arguments jl2| not involving any wire-specific 
information other than their aspect ratio {(/d) can be 
used to show that the latter quantity scales linearly with 
Vf{i/d), Vf being the volume fraction of the network, 
whereas («£;) scales with Vf{£/d)'^. The universality in 
these scaling laws means that films composed of different 
values of V/, £ and d may have similar connectivities. 

For the sake of generality and simplicity we choose 
to represent individual wires by 1-dimcnsional atomic 
chains within the tight-binding model with M-fold de- 
generacy in the orbital degrees of freedom. In this way, 
the Hamiltonian hj associated with a single wire labelled 
J corresponds to a linear chain that carries M quanta of 
conductance, i.e., hj = n' /i 1*^' J' A*!: where 
the state \n,j,fi) represents an atomic orbital fi localized 
at an atom numbered n within wire j and the sum over 
n and n' is for nearest neighbors only. The sum over 
/i ranges from 1 to M, which guarantees its M-fold de- 
generacy, and the parameter t defines the bandwidth of 
the individual wire. As we shall see, other Hamiltonians 
could have been chosen without major changes to our 
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conclusions. 

The Hamiltonian H associated with the network is not 
a mere sum of the individual hj over the index j but it is 
dependent on how they are connected to each other. To 
assemble the Hamiltonian of the full network we assume a 
number N of finite-sized atomic chains, each one of which 
containing 1.6 x 10"* atomic sites 16|. Bearing in mind 
that these chains ultimately represent NW of length I 
and diameter d, we can easily obtain the volume fraction 
Vf and the corresponding values for {a) and {ue)- The 
network is then created as follows. Chains are connected 
randomly so that there is a total of N{a)/2 connections. 
The intra-chain atomic sites that make the connections 
are also randomly chosen. The connection is introduced 
by a hopping probability 7 inserted in the Hamiltonian 
connecting two different chains, that is, by a potential 
V — \n,j,fi)'-f{n',j',^\. Attention is given to the fact 
that j and j' must necessarily be different to represent 
distinct wires. Notice that this contribution is diagonal in 
H indicating that the inter- wire coupling does not mix the 
M conducting channels carried by each wire. Regarding 
the role played by the electrodes, they act as charge reser- 
voirs and can be mimicked by using semi- infinite chains. 
Since (as) gives the average number of contacts to the 
electrodes (per unit area), we must include the appro- 
priate number of semi-infinite chains Ne — 2{aE) into 
the network. In summary, our model consists of a net- 
work composed of N finite chains randomly connected 
to each other plus Ne semi-infinite ones, each of which 
connected to a single finite-sized counterpart. The semi- 
infinite chains can be sub-grouped into two sets repre- 
senting contacts to the left and right electrodes. Further- 
more, no finite-length chains are allowed to be connected 
to more than a single semi-infinite one. This requirement 
is necessary to avoid short circuits. 

Having defined the network Hamiltonian, we can now 
use the Kubo formalism, to calculate the zero-bias con- 
ductance r across the film [3 [HI . Absolutely equiv- 
alent to the Landauer formalism, the Kubo formula for 
the conductance involves only a few Green-function (GF) 
matrix elements evaluated at the Fermi energy, which can 
be easily obtained by efficient inversion techniques [20j . 
One of the advantages of using this formalism is that any 
changes in the Hamiltonian are automatically accounted 
for in the corresponding GF. As we shall see, replacing 
the linear-chain Hamiltonian with another that describes 
the detailed atomic structure of NT becomes a straight- 
forward task, albeit more computationally demanding. 
The coupling between wires is regulated by the param- 
eter 7. Since we are interested in the upper bound for 
the conductivity, we can increase 7 until the conductivity 
saturates. At that point we can say, by inspection, that 
the conductance across any two wires connected by this 
value of 7 matches the intrinsic conductance of the wires, 
i.e, MFq, where Fq is the quantum of conductance. In- 
evitably, there will be some degree of reduction on the 
conductance when two NW are connected, particularly 
in the case of finite sized wires but the best way to min- 
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FIG. 1: Upper bound for the conductivity CTu (in units of S/m) 
as a function of the quantity x — Vft/(f (in units of m^^). 
Line with square symbols is for an array of linear chains, each 
one of which capable of carrying 2Tq of conductance. Circles 
correspond to results for disordered arrays of finite-sized NT 
in which the details of their electronic structure have been 
fully taken into consideration. In this case, 7 = 70, where 
70 is the tight-binding electronic hopping between nearest- 
neighbor carbon atoms of the NT. Error bars account for the 
standard error in the configurational average. 



imize this reduction is to consider very large wires. In 
the case of atomic chains, 7 = Hs an excellent choice for 
the coupling parameter since it reduces the conductance 
between two chains only by a very small fraction. 

With all the parameters defined, the upper bound for 
the network conductance can now be calculated for dif- 
ferent volume fractions and aspect ratios. Each calcu- 
lation involves a large configurational average in order 
to achieve statistical significance. Rather than the calcu- 
lated conductance values, we present in Fig. 1 the respec- 
tive conductivities plotted as a function of x = Vft/d?. 
With NT in mind, the number of conducting channels 
is M = 2. A wide range of lengths, diameters and vol- 
ume fractions have been considered but, because of the 
universality behavior of the connectivities, they all fall 
onto a common straight line (square symbols) described 
by cr„ = 13m X, where Pm = M x 4.25 x 10"^ S. This 
indicates that the dominant factor in determining the 
conductivity of the network is clearly the number of con- 
nections with the electrodes, which also scales linearly 
with X. 

The expression above provides the upper bound for 
the conductivity of disordered networks in the case of 
wires capable of carrying M quanta of conductance. It 
is instructive to test the expression for typical NT val- 
ues, namely £ = l^m and d — 1.2 nm. For Vj = 30%, 
the predicted upper bound would be au — 1-8 x lO^S/m 
if NT of these dimensions could be fully dispersed to 
form the network. NT are however known to bundle 
together, which means that in reality wire diameters 
are considerably larger. On the other hand, larger- 
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diameter bundles have more NT on the surface lead- 
ing to more current-carrying channels per wire. Taking 
all this into consideration, we can compare our expres- 
sion with the highest-conductivity case reported so far 
{Vf = 30%, £ = 5fim,d = 20 nm) [1 111 IH- Our 
prediction of CTu = 9 x 10^ S/m is only one order of mag- 
nitude superior to the measured value of cr = 6 x 10^ S/m. 
Bearing in mind that the upper bound here obtained as- 
sumes a number of ideal conditions that are experimen- 
tally unavoidable, this might be a clear indication that 
we are approaching a saturation point in the conductivity 
of NT-network films. 

No qualitative change is observed when the linear- 
chain Hamiltonians hj describing the individual NW is 
replaced with another that accounts for the precise elec- 
tronic structure of NT. In this case, a straight line with 
circular symbols depicts the conductivity of a network of 
realistic NT. The difference in slope and the fact that 
individual values are somewhat more staggered than the 
linear-chain results is easily explained by finite-size ef- 
fects. In the case of hj describing real NT, the number 
of atoms required to generate a wire of similar length £ 
is considerably larger. This introduces undesirable fluc- 
tuations in the conductance of individual NT, which on 
average lowers the overall conductivity. By making the 
wires longer, the fluctuations are reduced and the re- 
sults for NT tend to approach those for the linear chain. 
While we cannot increase the NT size without paying 
hefty computational penalties, we can reduce the size of 



the liner chains for the sake of comparison. In doing so, 
the linear behavior observed for linear-chain Hamiltoni- 
ans coincides with that for NTs. This suggests that the 
line with square symbols is therefore the most represen- 
tative for the upper bound for the conductivity of a NT 
network. 

Finally, although our focus has been on disordered net- 
works comprised of NT, we can extend our results to deal 
with other wires. This could represent the case of net- 
works made of other conducting materials, such as noble- 
metal wires, for instance. In this case the number of con- 
ducting channels M depends linearly on the wire diame- 
ter and the overall conductivity of the network is likely 
to scale inversely with d, with a proportionality constant 
that depends on the specifics of the wire in question. If 
the conductivity of NT-network films is approaching its 
saturation point, it is likely that wires other than NT 
may occupy the post of ideal components for disordered- 
network films. 

In summary, we have calculated the upper bound for 
the conductivity of NT-network films by assuming ideal 
contacts between tubes. Our results may be used to indi- 
cate how much room there is to lower the junction resis- 
tance within a film. More importantly, when compared 
with the highest measurements reported for NT-network 
films, the upper bound presented here points to a situa- 
tion in which the conductivity is approaching its limiting 
value. 
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